Conventional mechanisms for "exotic" superconductivity 
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We consider the pairing state due to the usual BCS mechanism in substances of cubic and hexag- 
onal symmetry where the Fermi surface forms pockets around several points of high symmetry. We 
find that the symmetry imposed on the multiple pocket positions could give rise to a multidimen- 
sional nontrivial superconducting order parameter. The time reversal symmetry in the pairing state 
is broken. We suggest several candidate substances where such ordering may appear, and discuss 
means by which such a phase may be identified. 



Most conventional superconductors are described very 
well by the BCS theory [y. The electron-phonon inter- 
action mediates an attraction between electrons that is 
stronger than the Coulomb repulsion. This gives rise to 
the Cooper instability of the normal state leading to the 
appearance of a condensate of pairs. The order param- 
eter (the anomalous function JF [Q) in this case belongs 
to the "s-wave" type, i.e., it is invariant with respect 
to the transformations oi G (S) R, where G is the crystal 
point group and R is time reversal operation. As a re- 
sult the quasiparticle spectrum has a gap, which leads to 
well-known experimental consequences. A variety of ma- 
terials: ^He |], UBei3 g, UPtg §, high T^ materials [|, 
and Sr2Ru04 m, have been discovered that potentially 
break the G (E) R symmetry of the normal state. A well 
known example is ^-phase of "^He g which is not rota- 
tionally or time reversal invariant (note that the B phase 
of ■^He is both rotation and time reversal invariant). Such 
non s-wave superconductors are usually expected to have 
a gapless excitation spectrum and arise when the inter- 
action itself depends upon the superconducting ground 
state (in the case of '^He the BCS ground state is the B 
phase and the spin fluctuation feedback effect is required 
to stabilize the A phase [g[). All the possible symmetry 
classes of the superconducting state in crystalline mate- 
rials were enumerated in Ref. [^ (for a review see Ref. 

00)- 

We show below that exotic superconductivity can be 
a much more common phenomenon and does not require 
unusual mechanisms. The electron-phonon and Coulomb 
interactions are enough to give rise to an multidimen- 
sional order parameter which would have lower symmetry 
than the normal state - including the breaking of time- 
reversal invariance. The effects we consider are possible 
in metals with several pockets which are centered at or 
around some symmetry points of the Brillouin Zone (BZ). 
A BCS approximation generalized to the the nrulti-band 
case (see, e.g., pj]) will be used. The new point here 
is that since the form of the interaction parameters de- 
scribing the two electron scattering on and between the 
different pockets of the Fermi surface (FS) is fixed by 
symmetry, the resulting superconducting state need not 



be s-wave. Below we consider three cases in detail: a) 
three FS pockets centered about the X-points of a simple 
cubic lattice; b) three FS pockets at the M-points of the 
hexagonal lattice; c) four FS pockets at the L-points in 
the face-centered cubic lattice. A complete analysis of 
all other high symmetry points is possible, and will be 
published elsewhere Q. 

We emphasize that this FS structure is not unusual. 
Indeed, superconductivity with pockets as in case (a) 
and Tc ~ O.IK is found in LaBe ||lj]. Another example 
is given by the superconducting semiconductors such as 
PbTe, SnTe, or SrTiOs M^. Many materials exist where 
such FS sheets co-exist with other non-symmetry related 
FS sheets and some of these materials have anomalous su- 
perconducting properties. One example is CeCo2 |16| , p7| . 
Fig. |l| shows some of the FS sheets of CeCo2 (Tc = 1.6 
K). We will return to this later O. 




FIG. 1. Some FS sheets of CeCoa (from Ref. fTTl). 



We will use the generalized Ginsburg-Landau (GL) 
functional to identify possible nontrivial superconduct- 



ing phases. The Hamiltonian for several separate pieces 
of the FS can be written in the following form: 
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where u and a' are spin indices, \ap(c\) includes the in- 
teraction for scattering two electrons from the pocket a 
into pocket (3 which is due to both Coulomb and electron- 
phonon terms. Introducing the anomalous Green's func- 
tion Ta(x — x') for each FS sheet a, the corresponding 
Gor'kov equations ||l^ can be used to obtain the follow- 
ing solution at finite temperatures for the case of singlet 
pairing: 
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Eq.(2) being expanded in |Aq| up to the third order, 
Eq.(3) becomes a variation of the GL functional with re- 
spect to the vector order parameter Ac(p) |l^,|2^. For 
simplicity we assume that each Ac(p) is constant along 
the corresponding FS, and to fourth order in Aq, we can 
write 
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where A~^ is the matrix inverse to the interaction Aq.^, 
ijJd is the cutoff (Debye) frequency. 

We now analyze three different cases for multiple FS 
sheets. 

(a) Three X points in a cubic lattice. 
The interaction matrix A for three X points takes the 
following general form: 



Aq/3 = XSa[3 + fJ-il — Sail)- 



(5) 



Here A is the interaction on the same pocket, /i couples 
any two different pockets. Gonsider first the linearized 
gap equation Eqs.(g) and ^ to determine T^- 
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The three Aq transform among each other at cubic sym- 
metry transformations forming a 3D reducible represen- 
tation of the cubic group Oh, which is split into a ID 
Aig and a 2D Eg irreducible representation. These two 
representations correspond to different order parameters 
with two critical temperatures: 
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(the terms in the exponents must be negative for the 
Cooper effect to take place) . The basis wave function for 
ID identical representation is 



l^ (Ai + A2 -H A3)/%/3 



(9) 



and the basis wave functions for the 2D representation 
can be chosen as 
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where e = exp(27ri/3). From Eq.(g) if A — /x < and 
/i > then superconductivity will belong to the nontriv- 
ial 2D Eg representation, i.e., if the interaction between 
two different FS pockets is dominated by Coulomb repul- 
sion. Let us consider the latter case in detail. Rewriting 
the Landau functional Eq.(y) in terms of I, rji and 772, we 
obtain, for temperatures T near T^.e'- 
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where F^ ' is the fourth order term in the GL functional With Tc^a < Tc^e the superconducting instability will 
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which may admix the \D representation I: 
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correspond to the 2D representation Eg of the cubic 



point group. The fourth order coefficients in Eq.([ll|) in- 
dicate that the class 0{D2) is the most preferable en- 
ergetically p. This class corresponds to a phase with 
^72 = 0, rji y^ in Eq.(0). The symmetry properties 
of this class are known K. Time-reversal symmetry is 
broken and allows for antifcrromagnetic domains and for 
fractional vortices to appear (see, e.g., [Q). In princi- 
ple point nodes should appear where the FS intersects 
the cube diagonals. This would lead to the electronic 
contribution in the T^ behavior for the heat capacity at 
low temperatures. In our case, however, there is no FS 
along the diagonals of the cube and the low temperature 
thermodynamic properties will be determined by the gap 
of the same magnitude for all three FS sheets. In the 
presence of another FS, for example, at the F-point, the 
nontrivial order 0{D2) will be induced on it. In this case 
the point nodes will exist and power laws in thermody- 
namic properties due to the superconductivity should be 
seen experimentally. The anisotropy of the upper critical 
field (-ffc2) 12^ near T^ for this class also requires that 
(at least) two vortex lattice phases (with a second order 
transition between them) exists when the magnetic field 
is applied along the (1,1,0) and equivalent directions. 
Note from Eq.(|lJ) that the terms linear in / identically 
disappear for this class, i.e., there will be no admixture 
of the s-wave component. 

(b) Three X points in a hexagonal lattice. 

Calculations in this case are the same as in the pre- 
vious case, i.e., Eqs.([5|)-(p^) apply. We only have to 
specify the symmetry properties and the superconduct- 
ing class in this case since the symmetry of the lat- 
tice is different. The three-dimensional representation 
is split by the hexagonal group D^h (^ R into lD{Aig) 
and 2D(i?2g). Note that the basis functions for £'23 can 
be once again chosen as given by Eq. (|lO|) . The phase 
with 7]i 7^ 0, ?72 = has the lowest free energy and 
in this case corresponds to the nontrivial class Dq{C2). 
Time-reversal symmetry for this class is broken M , and 
ferromagnetism is allowed. Point nodes (at two points 
of intersection of an additional FS at the F-point with 
the six-fold axis) can be seen in thermodynamic prop- 
erties but again these nodes are only present if such a 
FS exists. The upper critical field is isotropic near Tc 
for this superconductivity representation. Nevertheless, 
it can be shown that there will also exist (at least) two 
distinct vortex lattice phases with a second order transi- 
tion between them for the magnetic field applied in the 
basal plane | p^ . 

(c) Four L-points in the fee lattice. 

The interaction and the linearized gap equation for the 



four L points again take the form Eqs.(|2-|6|). This time 
the 4D representation A^ is split into the ID(Aig) and 
3D(i^2g) irreducible representations. 
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FIG. 2. Regions of existence of different superconducting 
phases on tlie basis of the three-dimensional representations 
of the cubic group (from Ref. pi). A point on the boundary 
of two phases corresponds to our BCS solution. 

The critical temperatures for the two representations 
are now given by 
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The basis functions for the ID and 3D representations 
are 



/ = (Ai + A2 + A3 -H A4)/2 (ID) 



(15) 



and 



r,^ = (Ai - A2 - A3 -1- A4)/2 

r?j, = (Ai + A2 - A3 - A4)/2 {3D) (16) 

T], = (Ai-A2 + A3-A4)/2. 

The nontrivial 3D representation is stable if A — /i < 
and fi > 0, i.e. if the interaction is attractive for each 
pocket alone, while it is repulsive between two different 
pockets. As above, we can expand Eq.(Q) in terms of I 
and 77. Dropping the ID identical representation, we get: 
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The GL coefficients in Eq.(|l7|) piaces tlie system riglit 
on tlie boundary of two phases, superconducting classes 
D^a\D2) <» R and DiiE) (see Fig.|). This degeneracy is 
an artifact of the BCS theory, it is not lifted by higher 
order terms in the GL functional. The presence of a FS 
at the r point lifts this degeneracy. As a result, the mag- 
netic superconducting class D4^[E) is likely to appear. 
This class also allows ferromagnetism. Note that it has a 
line of nodes (at the intersection of the FS with the hor- 
izontal plane of symmetry) i.e., C ocT^ at low tempera- 
tures 1^. There are also point nodes at the intersection 
of the FS with the four-fold symmetry axis. This repre- 
sentation also exhibits an anisotropy of Hc2 near Tc. In 
principle multiple vortex lattice phases can also exist for 
this class but they are not required by symmetry as they 
are for the E representations discussed above. 

In the above cases, other than the standard isotropic 
order parameter, only multidimensional order parame- 
ters appeared. This leads to the possibility of domain 
walls between different equivalent superconducting states 
and as a possible consequence the existence of inhomo- 
geneous magnetic order |^ (see also [ p^|Jll| ). Also in 
all the above cases the resulting superconducting states 
had gaps of equal magnitude on each of the FS sheets. 
In such a case a Hebel-Slichter peak in 1/Ti measure- 
ments may be present. The non-trivial representations 
were stable when the pair interaction between the differ- 
ent sheets was repulsive (independent of the intra-sheet 
interaction). In many materials the Coulomb repulsion 
can be comparable to the attraction between electrons 
due to electron-phonon interactions. This is illustrated 
through the reduction of isotope effect due to Coulomb 
repulsion. It is well known that for a number of metals 
Tc oc M~", where not only a 7^ 0.5 but it may even have 
the opposite sign, (a ~ —5 for a — C/ |g^). The sensi- 
tivity of these exotic superconducting phases to impuri- 
ties needs a more detailed analysis. However, provided 
the inter-pocket defect scattering amplitudes are much 
smaller than the intra-pocket amplitudes, these phases 
will survive the presence of a considerable amount of de- 
fects (due to the ordinary BCS pairing on each sheet). 
This will be studied in more detail in |13[ . 

In summary, we have shown that exotic superconduc- 
tivity can appear merely as a competition of the phonon 
and Coulomb interactions if the FS consists of several 
pockets located at some symmetry points. Time-reversal 
symmetry is broken for the nontrivial order, meaning 
that the superconducting transition should be accom- 
panied by some kind of magnetic order. The simplest 
methods to identify exotic order parameters are, apart 
from the phase-sensitive measurements, the power law 
dependence of the heat capacity (due to the nodes), mea- 
surements of the upper critical field anisotropy Hc2 {&) at 
Tc UM , or the observation of transitions between differ- 
ent vortex lattice phases. Note that if the FS pockets 



are fully isolated then the nodes are absent since the or- 
der parameter is then constant on each FS pocket and 
changes phase as one moves from one pocket to another. 
Nodes could appear, however, if there are "necks" con- 
necting different sheets mm or if superconductivity is in- 
duced on a FS centered, for example, around the F-point. 
The upper critical field anisotropy near Tc does not work 
as a test of nontrivial order in the hexagonal group pi| . 
The magnetic order, on the other hand, can be observed 
in /iSR measurements or magnetization measurements in 
small enough samples (where the dimensions are on the 
order of the penetration depth) . A general classification 
of all cases of nontrivial superconductivity of the type 
considered above is possible. These FS sheets are not 
always centered on the BZ boundary, as, for example, 
in some doped semiconductors ||l5[ and CeCo2 [|7|- We 
postpone a detailed analysis of the various possibilities 
to a future work fl^ j. 
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